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FOR BASE-POSITIONS HAVING' MAJOR/MINOR 

BASE CALL DISCREPANCIES CALCULATE: 
SCAN POSITION DIFFERENCE, PEAK HEIGHT 
RATIO, PEAK AREA RATIO. 
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Table 1 

System components and tunable parameters 



Program 


runciion peirorrneu 


BlastParse.pl 


Parse BLAST output 


mark substitution heteroplasmy.pl 


A A _. _l * i-i in -I l-i till §^\.t l/\ r\. <^ f** f\ ^^ollO 

Mark ambiguous oase cans 


extract_SE_consensus.pl 


Extract a consensus sequence from a 
phrap poiypnrea run, replacing iow quality 
bases with '?', and polyphred rank 1 

baSBS Wltn IN . 


compute_coverage.pl 


Calculates overall, forward strand, and 
reverse strano coverage toi dn ua&c& m 
the consensus sequence 


count_hv1_deletes.pl 


Count deletions in HV1 relative to the 
rOKb 


border_index.pl 


Compute start and end positions of HV1 
and nvz regions. 


fix_mitotype_reporting_range.pl 


Output a mtDNA profile based on input 
variants list and reportable range. 


find_bad_traces_from_blast_report.pl 


Reject reads that do not align 
appropriately to me ror\o 
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Adds template name, template type, and 
orimer tvoe to phred output files. 


seq2delta_vs 


Align mtDNA profile to rCRS and report 
variants according to nomenclature. 
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Mask mtDNA consensus sequence based 
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coverage 
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Count insertions in HV1 relative to the 
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Extract the HV1 portion and HV2 portion 
of the mtDNA consensus sequence, 
based on the computed HV1 and HV2 
reqions. 
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Invert an X masked sequence output by 
cross match 
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Table 2 

Assessment 1 Assessment 2 



Total True Positives: NantnmateH = Nmanuai 


16 


10 


Total True Negatives: Automated = Amanuai 


77,358 


110,354 


Total False Positives: N au tomated = A manu ai 


95 


807 


Total False Negatives: A au tomated = N ma nuai 


6 a 


14 


Total Incorrect 


0 


0 


Sensitivity: TP/(TP+FN) 


72.73% 


41.67% 


Specificity: TN/(TN+FP) 


99.88% 


99.27% 


Positive Predictive Value: TP/(TP+FP) 


14.41% 


1.22% 


Negative Predictive Value: TN/(TN+FN) 


99.99% 


99.99% 
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